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Recent sum rule analyses on the <VV-AA> two-point correlator have led to significant discrepancies in the 
values found for the OPE condensates, most dramatically in the dimension eight condensate and to a lesser 
extent in the dimension six one £Q. Precise knowledge of these condensates is of relevance in kaon decays [21 and 
therefore it seems mandatory to assess the actual impact of what is commonly neglected in spectral sum rules, 
most prominently the issue of duality violations. We will explicitly compute them in a toy model and show that 
they are a priori non-negligible. 



1. INTRODUCTION 

In the absence of a solution to QCD, extracting 
information about QCD Green functions is a chal- 
lenging task. Experimentally, only their imaginary 
parts, i.e., the spectral functions 



PL R {t) = -lmU LB (t) 



(1) 



are accessible. A common strategy is to use 
Cauchy's integral theorem smartly, so that spectral 
functions can be related to the full Green functions 
through 
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and replace HL R (q 2 ) by its Operator Product Ex- 
pansion, defined as (Q 2 = —q 2 ) 
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However, since the OPE is not convergent over the 
whole <j 2 -plane (most blatantly over the physical 
axis) the previous replacement generates an error 
2?[™l(s ) so that Cauchy's theorem now reads 

f S dtt n p LR {t) = ~l dq 2 q 2n Ii° L P R E {q 2 )+V^{ So ). 
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From now on, we will denote each term in the pre- 
vious equation as 

M n (so)=A n {so)+V^( So ). (5) 

The moments M„(so) are extracted from exper- 
iment, whereas the OPE integral can be done 
straightforwardly to yield 
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The so) functions account for the difference 
between nifl(g 2 ) a.ndUl R E (q 2 ), i.e., they measure 
the amount of duality violations [3]. It is worth 
stressing that, even though present-day analyses 
do not implement them, as a matter of principle 
f ' n '( s o) 7^ 0. In the next section we will attempt 
to characterize such terms. 

2. WHAT WE EXPECT ON DUALITY VI- 
OLATIONS 

Our discussion will rest upon very mild assump- 
tions, namely that the OPE is an asymptotic ex- 
pansion and that large- N c QCD 0j is a good ap- 
proximation to real QCD. In the deep Euclidean 
regime, a Green function can be approximated by 
its OPE up to exponentially suppressed terms (the 
typical uncertainty in asymptotic expansions) 
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A being a typical scale. To be conservative, we 
consider the OPE to be valid all over the left half 
complex q 2 -plane. We can safely extend Q over 
the whole q 2 -plane (\q 2 \ 3> 0) but on the physical 
axis, where singularities are located. This means 
that over the right half g 2 -plane, except for the 
analytically-continued OPE terms there can only 
appear exponential pieces, which wc denote under 
A(q 2 ), and 



\l LR {q 2 )^U% E {q 2 )+A{q 2 )+0[e 
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This extra piece A(g 2 ) has to depend somehow on 
the number of colors N c . Recall that in the strict 
large- N c limit resonances are infinitely narrow and 
nothing like an exponential damping is expected 
over the physical axis. It turns out [H] that A(q 2 ) 
can be cast in the form 



A(q 2 ) ~ exp 



H(q 2 



(9) 



where !F{4>, N c ) is a function such that in the 
Nc — ► oo limit the exponential behaviour no 
longer holds along the Minkowski axis (0 = 0, 2tt). 
Regarding 7i(q 2 ), not much can be said in general. 
However, bearing in mind that Ulr is the differ- 
ence between vector and axial currents, at least its 
imaginary part has to be of an oscillatory type. 
Comparing JHJ) and one can see that A(g 2 ) and 
f ( s o) are intimately related. Actually, 

pW(So) = ~i / ^ A ^+° 
Re q 2 > 

(10) 

Use of Cauchy residue theorem bearing in mind the 
exponential behaviour of © allows one to eventu- 
ally express the previous equation as 2 



dt t n - Im A(t+ie) + O (e~^\ 

(11) 



i.e., only in terms of the imaginary part of @. 

3. A TOY MODEL 

Lacking knowledge of QCD Green functions, in 
order to proceed further one has to resort to mod- 

2 See [S] for a detailed derivation. 



els, where analytic results are available. Our aim is 
two- fold: 1) we will be able to test different spectral 
sum rule techniques and 2) do an explicit calcula- 
tion of the duality violating pieces. 

We will consider the following ansdtze for the 
spectral functions 



Pv 



PA 



2F 2 S(t - m 2 ) + 2 FvW ~ M vH) 

oo 

2F 2 S(t) + 2j2F 2 A S(t^M 2 A (n)) 7 (12) 



n=0 



where both axial and vectorial towers display a 
Reggc-like behaviour 
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For simplicity, we take a single spacing A 2 and 
decay constant F 2 . So far, our model displays a 
spectrum of infinitely narrow vector and axial res- 
onances. We want to provide them with a non-zero 
width while preserving the analytic properties of 
nii?(q 2 ). A simple choice |J] is to shift q 2 to 
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With the previous replacement, widths have the 
expected A c -scaling 
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and the full Green function reads 
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Free parameters are fixed by demanding the short- 
distance behaviour of QCD to be matched, namely 
the parton model coefficient for both vector and 
axial channels and the two Weinberg sum rules. 
From l|l()|> it is straightforward to obtain the OPE 
coefficients. To leading order in a/N c they read 
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Bernoulli polynomials display (asymptotically) a 
factorial growth. This signals at an asymptotic 
OPE for our model, as it is believed to happen in 
real QCD. 

The game we want to play now is the following: 
we will take our toy model as the real world and test 
the different existing approaches, mainly FESR and 
pFESR |B] upon it. FESR rely on the assumptions 
that Aq,i = and that there exist duality points 
Sq which lie rather close for neighbouring moments, 
i.e., 

(18) 

On the other hand, the pFESR approach relies on 
a slightly modified version of J2J, namely 

dtw{t) PLR {t) = --^ I dq 2 w(q 2 )U LR {q 2 ), 

10 2m J\q 2 \=s 

(19) 

where w(t) are polynomials which, because they 
vanish for t = ,s , arc expected to minimize dual- 
ity violations. In pFESR, one has a combination of 
moments that has to match the corresponding com- 
bination of OPE coefficients. Usually, such OPE 
coefficients are fixed by fitting both curves over a 
window in Sq. 

However, once applied to our model, both FESR 
and pFESR yield predictions with huge errors when 
compared to the true OPE coefficients. We defi- 
nitely have to do better and assess the impact of 
the neglected terms. If we take the imaginary part 
of ©, 



Eq. (|22[1 is nothing but a particular example of 
what we already stated in @, where Im Ti.(q 2 ) 
is, to a very good approximation, of a sinusoidal 
type. Eq. (|21[1 is the analytical continuation of the 
OPE into the Minkowski region. Therefore, the ab- 
sence oft -1 and t~ 2 powers is a direct consequence 
of the Weinberg sum rules we enforced earlier on. 
Those power corrections are closely related to the 
second line of © and have a (in principle) non- 
negligible impact, presumably already at dimen- 
sion eight, and surely for higher order condensates, 
where they have positive powers of sq out front. 
Both duality violations and power corrections I© 
are generic and nonetheless missing in all existing 
analyses pQ. Any strategy for improvement should 
take them both into account. 

One possibility J5] is to take the first two mo- 
ments in whose explicit form is 
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Mi (so) = 
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with Im A(q 2 ) as given by (|22|l 3 and determine the 
parameters k, 7, a and (3 through a fit. In the 
above equations, terms coming from the OPE are 
strongly suppressed: first, Weinberg sum rules en- 
sure that C2 and C4 vanish, and it seems reasonable 
to neglect (at first) the remaining 62fe-terms due to 
the suppressing powers of s . Then a prediction for 
A2 3 is possible through 



Imn ifi (t) = Imn^(t) + ImA(t), 



(20) 



and use our model to explicitly compute the differ- 
ent pieces we end up with [S] 
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M 2 (s )+ dtrimA(t) = 
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for a given sq. Notice that this strategy shares 
some similarities with FESR, but we need no longer 
assume the same duality points for different mo- 
ments: we can determine the OPE coefficients, e.g., 
through a fit in a window of So- In order to make 

3 The use of a sinusoidal function instead of the (unknown) 
Im Ti.(q 2 ) is a justified approximation as long as one does 
the fit over not too large a window. 



a comparison, we list on the table below the re- 
sults found in our model, where one can see that 
the duality violation approach (DVA) we propose 
outdoes conventional FESR and pFESR. 





A 6 {GeV*) 


A 8 (GeV») 


FESR 


-2.0- 1CT 3 


-1.6- 10~ 3 


pFESR 


-3.8 • 10" 3 


+6.5 • 1CT 3 


DVA 


-2.51 • KT 3 


+3.29- 10~ 3 


OPE 


2.81 10 3 


+3.44 10 3 



4. CONCLUSIONS 

Based on rather mild assumptions, namely that 
the OPE is an asymptotic expansion and that the 
1/N C expansion is a good approximation to QCD, 
we would like to propose a general expression for 
the duality violating pieces in Hlr(o 2 )- With the 
help of a toy model inspired by large- ./V c we have 
been able to compute duality violations explicitly. 
Then, using the model as if it were the real world, 
we have tested FESR and pFESR upon it. Both 
approaches yielded predictions for the condensates 
with large errors, showing that duality violations 
might have a bigger impact than commonly ex- 
pected. We propose an analysis that does include 
duality violations and in which errors reduce down 
to typically 10-15%, which is a clear improvement 
over FESR and pFESR. Interestingly, the method 
can also be applied to QCD, since it relies on 
very generic assumptions. Finally, we would like 
to emphasize that duality violations are not the 
only source of uncertainty in present-day analyses: 
one should not forget about potentially harmful so 
power corrections (see second line of ©), whose 
impact might be non-negligible already at the di- 
mension eight condensate. 
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